Incoherent dynamics of vibrating single-molecule transistors. 
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We study the tunneling conductance of nano-scale quantum "shuttles" in connection with a recent 
experiment [H. Park et al, Nature, 407, 57 (2000)] in which a vibrating Ceo molecule was apparently 
functioning as the island of a single electron transistor (SET). While our calculation starts from the 
same model of previous work [D. Boese and H. Schoeller, Europhys. Lett. 54, 668 (2001)] we 
obtain quantitatively different dynamics. Calculated I-V curves exhibit most features present in 
experimental data with a physically reasonable parameter set, and point to a strong dependence 
of the oscillator's potential on the electrostatics of the island region. We propose that in a regime 
where the electric field due to the bias voltage itself affects island position, a "catastrophic" negative 
differential conductance (NDC) may be realized. This effect is directly attributable to the magnitude 
of overlap of final and initial quantum oscillator states, and as such represents experimental control 
over quantum transitions of the oscillator via the macroscopically controllable bias voltage. 

PACS numbers: 73.63-b,73.23Hk,72.20Dp 



INTRODUCTION 

Over the past decade, nanotechnology has advanced 
the ability to build systems in which chemical self- 
assembly defines the functional and structural units of 
nanoelectronic devices 0, IS IS 01 • Since elastic pa- 
rameters of organic compounds currently utilized in such 
work can be much " softer" than those of semiconductors 
or metals, one must carefully consider the important role 
mechanical degrees of freedom may play in charge trans- 
port devices |S|- It is also fundamentally interesting to 
consider how the physics develops as one proceeds from 
the nanoscalc all the way up to macroscopic electrome- 
chanical systems |9|, [lOj . 

A significant impetus in the field of chemically self as- 
sembled nanostructures is in the direction of Coulomb 
blockade devices, in which transport is due to the tunnel- 
ing of a single electron or Cooper pair through the device 
via island structures of small capacitance [S Ull ll2L fl3| . 
The quest for systems with ever smaller capacitance has 
driven workers to define smaller structures in such de- 
vices, with the result that current experimental work, 
such as that modeled here, involves the utilization of sin- 
gle large molecules as the island structure through which 
electronic transport takes place. 

The experiment of Ref. Q involves the fabrication of 
metallic gaps roughly 1 nm in separation |l4| . and ex- 
posure to a solution of Cqo molecules. When IV curves 
were measured, some of the gaps exposed to Cqo dis- 
played a robust 100 meV Coulomb blockade with charg- 
ing energy modulation via a third gate terminal (in this 
case the gate terminal was the underlying semiconduc- 
tor substrate). Also observed was a series of smaller 5 
me V- wide steps attributable to coupling of the electronic 
transport through the island to the quantized vibrations 
of the center of mass motion of the molecule in its Van der 
Waals binding potential nearest to one of the leads. The 



steps appear to have some remarkable secondary features, 
i.e., asymmetries about zero bias voltage, and unresolved 
structure in the intra-step regions of the IV characteris- 
tics. 

It was suggested in Ref. [l| that the number of vi- 
brational steps and their amplitudes are controlled by 
Franck-Condon physics. This idea was further developed 
in Ref. [S where Franck-Condon matrix elements were 
used to calculate transition rates between different vi- 
brational states and to solve the master equation for the 
I-V curves. Our work differs from that of Ref. in the fol- 
lowing respects: i.) While we work with the same model 
as Ref. y, we arrive at a different functional form for the 
transition matrix elements, with the result that our cur- 
rent voltage characteristics disagree quantitatively, ii.) 
Because of the unitary nature of the space-shift operators 
that must act on the electron tunneling perturbation, we 
stress the distinction of shuttle transport into unitary and 
non-unitary limits, depending on the degree to which en- 
ergy conservation allows the tunneling perturbation can 
fully connect the Hilbert spaces of states before and af- 
ter the island is charged, iii.) We also argue that a bias 
voltage dependence of the Franck-Condon factor should 
be considered, and that in fact such a dependence may 
explain the aforementioned secondary features in the IV 
characteristics of Ref. [l|. A stronger dependence is fur- 
ther shown to lead to a "catastrophic" negative differ- 
ential conductance effect which we will describe in this 
paper. The Chalmers group has modelled the system 
with classical, damped oscillations and incoherent elec- 
tronics 8], and with classical undamped oscillations and 
electronic coherence |2j . Our model is that of incoherent 
electronics and quantum mechanical oscillations that are 
strongly damped. As compared with the model of Ref. 8, 
the molecular island displays its vibrational character not 
by shuttling charge at a rate u , but in a more "spec- 
troscopic" manner, as a sort of fingerprint in the charge 
transport physics. Strong dissipation of the oscillator's 
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motion plays a key role in our incoherent dynamics, in 
contrast to Q, which assumes undamped classical mo- 
tion. 



THEORETICAL MODEL 

Fig. |QJ shows a depiction of the Cqq transistor sys- 
tem. The physics is assumed to be that of a 3 termi- 
nal quantum dot (including a capacitively coupled gate 
terminal) vibrating in a 1-d quantum harmonic oscilla- 
tor potential (we disregard vibrational states in other 
directions by assuming that their frequencies are much 
higher). We write the system Hamiltonian as follows: 



H = H, 



E 6 > 



Hi 



s=±l 



e 2 (h — n ) 2 
2C 

y(bath) 
int 



nEx 



H 



Bath 



(1) 



The vibrational portion of the island Hamiltonian is that 
of a simple harmonic oscillator in terms of the island 
center of mass momentum p and coordinate x 
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The right (s = 1) and left (s = —1) metallic leads are 
described by H s as normal Fermi liquids with constant 
Fermi-surface density of states pf, and chemical poten- 
tials controlled by the bias voltage. Charge transport 
between the island and leads is treated within the tun- 
neling Hamiltonian approach 
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where c\ 



and c\ a are the electron creation operators in 
the leads and island respectively. Although spin is con- 
served during all tunneling events we need to keep track 
of island electronic spin states, since tunneling events in 
the two spin channels are correlated. Here t is a constant, 
position independent tunneling matrix element, with the 
island centered between the metallic leads. Since the is- 
land is free to move, we can also expect an exponential 
dependence of tunneling on the center of mass coordi- 
nate; the decay of the conduction electron wave function 
in vacuum is given by the parameter A. 

Next we introduce the charging (Coulomb) energy of 
the island with capacitance C (fourth term) and its de- 
pendence on the island coordinate (fifth term), written 
in the form of charge-electric field coupling. The value 
of n a is controlled by the gate voltagejlg. In the present 
study we will assume that C is small, and n is very close 
to 1/2, so that at all relevant temperatures and bias volt- 
ages only two charge states, n = and n = 1, contribute 
to the dynamics. We denote the charging energy differ- 
ence Ec — (e 2 /2C)[(l — n ) 2 — n 2 ,]. Since the amplitude of 
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FIG. 1: Schematic of Ceo transistor system in terms of com- 
bined electronic-vibrational transitions for passage of current 
through the device. Dashed arrows indicate transitions in 
which quanta of oscillator energy hu) are supplied to the Cqo 
by the tunneling electron 



zero-point fluctuations, u — y/h/2muj , is rather small 
for the heavy island (it is estimated as ~ 4 pm for the 
Ceo molecule 0) we do not consider terms no n- linear in 
x. 

Finally, the last two terms in the Hamiltonian de- 
scribe an environment and it's dissipative coupling to 

the molecule via y(bath) ^ ^wo mechanisms may con- 
mi . . J 
tribute here. The vibrational state of the island cou- 
ples to the bulk and surface vibrations of leads via stan- 
dard displacement-dependent elastic energy. On another 
hand, island oscillations produce electrostatic potential 
fluctuations in the leads which excite electron- hole pairs. 
We consider the second mechanism as dominant because 
(i) it is Ohmic dissipative while the phonon coupling is 
not and vanishes for small energy transfer and (ii) 
it is mediated by the unscreened Coulomb interaction for 
the charged island. By writing the coupling term as 
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we arrive (using golden-rule approximation) at the bath- 
mediated transition rates between vibrational states of 
the island, \l >, 
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in the form characterized by a single parameter K cx 
(pfV) 2 . [In a more refined theory one may generalize to 
consider K depending on the island charge as well.] We 
do not see any obvious reason why this parameter may be 
anomalously small for leads with metallic concentration 
of conduction electrons. Thus transitions rates between 
molecular levels of order 10 _1 w o or 10~ 2 w o must be con- 
sidered as typical. We performed all our simulations for 
if = 0.1. 
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In what follows we consider the case of weak conduc- 
tance when tunneling rates given by r(°) = 2np F t 2 arc 
much smaller than KTiuj . This condition makes the the- 
ory very simple conceptually because one may ignore in- 
terference between subsequent tunneling events — on time 
scales larger than K/W the density matrix of the is- 
land reduces to its diagonal form in the oscillator eigen- 
states representation and one may discuss dynamics in 
terms of probabilities, P n ,a,i, of finding a system in a 
state \n,a,l >. The I-V curve is then obtained from 
the steady-state solution of the master equation. This 
rather standard framework is identical to the one used 
in Ref. All transition rates for the master equation 
may be calculated exactly for the model given above. In 
what follows we explain how this is done and analyze 
various limiting cases by solving equations numerically. 
Other work [3] considers the possibility that electronic 
coherence persists over timescales larger than ~ 
while assuming undamped classical motion of the oscilla- 
tor degree of freedom. By contrast, our master equation 
approach has strong dissipation of the oscillator motion 
as a built-in condition. In Q it appears that the stable 
limit cycles of oscillation are a result of a kinetic stability, 
and not energetic equilibrium of the oscillator with the 
environment. 

To calculate tunneling matrix elements we have to take 
care of final-state interaction effects originating from the 
nEx term, in other words we first have to diagonalizc 
the Hamiltonian in the absence of tunneling separately 
for n — and n = 1, and calculate matrix elements be- 
tween the eigenstates corresponding to n values differing 
by ±1. The resulting n = and n = 1 Hamiltonians 
are related to one another by a shift in space of the po- 
tential minimum. Technically, it is convenient to calcu- 
late matrix elements in the same basis, and account for 
the difference between initial and final Hamiltonians by 
constructing a unitary transformation relating the cor- 
responding eigenfunctions. This is the essence of the 
Franck-Condon principle. In our case we have to deal 
with the shifted-oscillator basis set which is a textbook 
problem. Oscillator states with n — 1 and n — are 
related by 
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The ground state energy levels are shifted by —j 2 ~huj . 

Since the oscillator Hamiltonian suddenly changes dur- 
ing a tunneling event, the above unitary transformation 



is all we need in order to calculate matrix elements be- 
cause it projects n = 1 states into the n = Hilbcrt 
Space. An effective form of the tunneling Hamiltonian 
(in the n — representation) is given now by 
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where dimensionless a = u Q /X measures the ratio be- 
tween the molecule zero-point displacement and electron 
localization length. The energy shift term —j 2 htu must 
be taken into account when writing the energy conserva- 
tion law for all transitions. Specifically, when an island 
charges/discharges, the tunneling electron will gain/lose 
an energy equal to the shift. The above considerations 
are readily generalized for any relevant values of charge 
states. 

So the picture we have is that of the combined 
electronic-vibrational system with the tunneling Hamil- 
tonian giving dynamics to both degrees of freedom. In 
the realistic limit of W ^> r™, we can assume incoherent 
dynamics and move to a master equation formulation of 
the the occupation probabilities of quantum states of the 
transistor: 
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By considering transitions through only one side of the 
island (e.g., the right hand side) we obtain the current as 
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We have two kinds of transitions in Eq. 1)10(1. the 
thermal rates from (0), which change only the vibra- 
tional state (n' = n), and the charge-transfer rates, 

r = r^+rw, with 
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(/, n',(j', l'\H T \i,n,a, I) 



xSiEf-Ei + en'^A'-en^i); {n'^n), (12) 

where n, a and n' , a' are the initial and final electronic 
state indices for the island, I, V are likewise the initial and 
final vibrational state indices, and i, f stand for the initial 
and final states of the leads which are considered to be in 
the state of thermal equilibrium (pf q ^ is the equilibrium 
density matrix of the leads). We also assume the leads 
to be symmetrically voltage biased, i.e. +V/2 on the left 
lead and — V/2 on the right one, to avoid unnecessary 
complications. For the same reason we do not consider 
the case when tunneling rates r^ ) are different for the left 
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and right leads. Explicit expressions for charge-transfer 
rates are: 



■ 0,0,/->l,<T,i' 

-.00 

l,<T,i'->O,0,i 

<5i; 



o|,7 (13) 
r(°) jo^j 2 [1 -/(«?)] (w) 

E c +'-^- + huj (l' -l -l 2 ) , (15) 



(where f(x) is the Fermi distribution function) with os- 
cillator matrix elements 



0$ = 0$ = e-^+^^-^VUW. 

min(l,l') . a'-m , \l-m 

(I - m)\(l' - m)\m\ 



(16) 



For all simulations we measure energies and voltages in 
units of 1 meV, and current in units of er^ ). 

In calculating IV characteristics from Eqns. (|l(Jfl and 
(|11[1 . we have to adopt a model for the bias voltage depen- 
dence of the electric field E at the position of the island, 
see Q. Phenomenologically, one might assume there is 
bias voltage independent component, Eg, which includes 
effects of the gate, charged impurities, image charges and 
screening, device geometry, etc., and bias dependent part 
proportional to V/d where d is the separation between 
the leads. Since the electric field for zero and non-zero 
V is calculated for the same charge and position of the 
island, we do not see any physical reason why the de- 
pendence of the field on V might be non-analytic. While 
this effect was implicit in the classical oscillator of 0, it 
was not considered in 0, . We argue below that it is 
precisely this term that may help to explain some of the 
anomalous features in the data of 0], and that negative 
differential conductance (NDC) may arise in systems in 
which this terms dominates over Eq. Indeed, this bias 
voltage dependence, along with gate voltage dependence 
of Eq, represent macroscopic "knobs" with which to con- 
trol the quantum mechanical transition matrix element 
between initial and final oscillator states. Thus we have 
to assume that coupling parameter 7 is a function of V, 
and write it as 



7; 



eEu Q 
huj n 



= c\ + (eV/hu; ) c 2 



(17) 



where c\ and C2 parameterize the relative contributions 
of bias voltage independent and bias voltage dependent 
electric fields. 



SIMULATION RESULTS 

Since for heavy islands we expect a <C 1 to be a typical 
situation, we start the analysis from the a = case first. 
When voltage dependence of 7 is ignored (as in Ref. 0), 
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FIG. 2: IV characteristics for c 2 =0,Ci = 0.0, 0.8, 1.5 with 
the voltage offset for each curve defined by the parameter 
E c — 10 meV (see text) 
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FIG. 3: IV characteristics for finite C2 = 0.005, c\ 
0, 0.8, 1.5, and E c = 10 meV. 



i.e. C2 = 0, we expect that in the limit V/hto >> j 2 all 
I-V curves must saturate to the ideal value Iq = 2eT^/3. 
The reason for that is in the unitary nature of Franck- 
Condon factors; for zero a we have 0\f] =< l'\e s \l >, 
and thus when all excitation channels are contributing to 

the current, we have J2i> \Om\ 2 = 1- Fig. © shows cal- 
culated IV characteristics in the strong dissipation limit 
for C2 = 0, ci = 0.0, 0.3, 0.8, 1.5. We have chosen (for 
all plots) hui = 5 meV, ksT = 0.15 meV, to correspond 
to the experimental data. 

There are several observations to be made. First, steps 
appear at roughly 10 meV intervals, corresponding to 
eV/2 = khuj , with k = 0,1,2,.... Note that it is only 
energy conservation in tandem with non-zero Franck- 
Condon transition matrix elements that gives multiple 
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FIG. 4: IV characteristics for a = 0.5, c2 = 0.05, 0.1, 0.25, 
and E c = 20 meV 



FIG. 5: Conductance resonance due to "catastrophic" nega- 
tive differential conductance effect for ci = 0.5, c2 = 0.5, and 
E c = 20 meV 



steps. Second, the current through the device saturates 
to the same limit for all values of c\. Also, since the 
action of the unitary operator is to project the initial os- 
cillator state (in the strong dissipation limit, this initial 
state is always the ground state of the oscillator) into the 
eigenbasis of a space-shifted Hamiltonian, we see that the 
larger the shift parameter 7, the more access the system 
has to transitions into states above the ground state of 
the shifted basis, resulting in more visible steps in the IV 
characteristic. This is in contrast to the results of Ref. 0, 
Fig. 3, where many steps are present when c\ = and 
the IV curves do not satisfy the unitary limit require- 
ment. The origin of discrepancy is hard to identify since 
the transition matrix elements were not given explicitly 
in Ref. |3j. 

We note, that for cy > 1 the structure of the curve de- 
pends on the value of the 2{E C — 'y 2 fnu )/fiuj parameter 
because different number of vibrational states are acces- 
sible for the electron to enter and to leave the island. 
These states appear at 2h~uj intervals and form two sepa- 
rate sets which do not coincide if 2(E c /frui — 7 2 ) is not 
integer. When one of the sets reaches its unitary limit 
before the other set starts conducting (as in the cases 
cl = 0.0, 0.3 in Fig. (J5J) we observe only the second set 
(thus curves with large offset show only one set of steps) . 
For small values of c\ conductance is dominated by just 
one or two states and the second set has no visible ef- 
fect on the curve ("fades away"). For the largest value 
of ci = 1.5 we see that neither set of steps (entering or 
leaving island) reaches a unitary limit before the other 
and hence both sets of steps are resolved. 

When voltage dependence is accounted for in 7, 
Eq. I|17|l . the I-V curves become asymmetric. The ob- 
vious effect of non-zero ci is that now the bias voltage 
affects the unitary transformation. In particular, we note 



negative differential conductance in the intra-step regions 
on the positive side of bias, while the intra-step regions 
at negative bias have positive differential conductance, 
see Fig. ©. This can be explained in the following way: 
in the intra-step region on the positive side of bias, as 
voltage increases, the shift parameter 7 increases, reduc- 
ing the overlap between the initial oscillator state and 
the final, shifted state (or states), and hence quantum 
mechanically constricting that particular I — ► I' channel. 
For negative bias voltages, we note that the constant elec- 
tric field parametrized by ci and the voltage dependent 
electric field parametrized by C2 work against one an- 
other, such that orthogonality between initial and final 
oscillator states is overall reduced. These I-V character- 
istics display what appears to be salient in the data of Q] , 
namely an asymmetry about zero voltage in differential 
conductance, particularly in the intra step regions. 

Fig. shows I-V curves for c 2 = 0.05, 0.1, 0.25, with 
ci = 0.5 and a = with charging energy E c — 20meV. 
Here we see what we refer to as "catastrophic" negative 
differential conductance. Although the tunneling oper- 
ator is still unitary, as ci increases, the energy-shift of 
the final oscillator basis grows strongly (as ~ —j 2 hui ), 
with the result that the initial state has its largest projec- 
tion onto states high in the final-state spectrum. As the 
bias voltage increases, the orthogonality effect is more 
and more important, until most of the (oscillator) over- 
lap is in a region where energy cannot be conserved, and 
so no channels for oscillator transitions remain, giving 
nearly zero current through the island. Finally, when 
\E C — (ci + (V/huj ) C2) 2 hu> \ gets larger than eV/2 the 
number of states available for the electron to tunnel off 
the island is reduced to zero. Fig. (J5J shows how this 
effect can produce a sharp conductance resonance. 

Fig. |JBJ shows the effect of a non-zero a. We have set 
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FIG. 6: IV characteristics with C2 = 0; Ci = 1.5 and a = 
0, 0.04, and E c = 20 meV 

a = 0, cl = 1.5, and C2 = for one of the plots, and 
a = 0.04, cl = 1.5, and c 2 = for the other, a = 0.04 
is a physically reasonable value which for the Cqq case 
corresponds to the localization length A ~ 1 A. Curves 
with finite a differ only slightly from the a — > case in 
the magnitude of current at various steps and asymme- 
try around V — 0. This is indicative of the manner in 
which a destroys the unitarity condition 1 2 = 1- 

Since we assume here that zero-charge island states are 
not shifted, the system conducts better when charged 
states are shifted towards lower voltage lead. We note 
that charge transfer physics discussed here radically dif- 
fers from the semiclassical limit of large- amplitude vi- 
brations presented in Refs.|U ,0 where electron tunneling 
was considered to be the fastest process at least for some 
values of x. In our case, finite a does not necessarily 
mean larger conductance at all voltages. Also, in the 
semiclassical limit of large-amplitude vibrations the an- 
swer crucially depends on dissipation — the vibration am- 
plitude is set by the balance between the energy gained 
from charge-transfer processes and energy lost to the en- 
vironment (in the limit of large c\ , c 2 and relatively large 
a). 



CONCLUSION 

Even at the level of the Hamiltonian one can easily 
see that the dynamics are controlled by many parame- 
ters. Although the master equation formally solves the 
problem for an arbitrary parameter set (as long as dissi- 
pation is strong enough), it is difficult to cover all pos- 
sibilites single paper, even with some of the simplifying 
assumptions used above. That said, we compare with 
the C 60 experiment and find that the IV characteristics 



can be more or less captured by the c\ ~ Lc2 ~ 0.05 , 
q ~ 0.1 parameter set. It does appear in [l| that there 
may be features (NDC in the intra-step regions) that re- 
flect the physics of the c 2 ^ case (see Fig. ©) when the 
space-shift in the oscillatory potential immediately after 
a tunneling event depends upon the bias voltage. It is 
intriguing to note that with c 2 ^ in our model, we have 
a situation in which the bias voltage, a macroscopically 
controllable variable, can have an observable influence on 
the quantum matrix elements which govern the dynamics 
of the interacting, quantum mechanical degrees of free- 
dom for the system. Future experiments on molecular- 
based transistors may further explore the existence of 
such effects, and point the way toward greater control 
over quantum transitions. The electronically coherent 
model 0], also seems able to reproduce primary features 
in the experimental IV curves, but come from otherwise 
qualitatively different behavior (an classically oscillating 
island with undamped, stable amplitude ~ 1 Angstrom, 
compared with our strongly damped, near zero point mo- 
tion (~picometers)). We find that the most interesting 
physics in our model is the (possibly catastrophic) neg- 
ative differential conductance effect arising from the ex- 
perimentally controlled Franck- Condon transition matrix 
elements. 
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